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© Surveying of boreholes. 

© Borehole surveying methods and apparatus for 
surveying the true longitudinal magnetic field within 
a substantially non-magnetic drill collar occupying 
the part of a borehole being surveyed, despite the 
collar being of insufficient length to provide longitudi- 
nal magnetic field measurements which are uncor- 
rupted by the longitudinal magnetic influences of 
adjacent magnetic drill string and bottom-hole as- 
sembly components. A plurality of longitudinal mag- 
netic field measurements are made by a static in- 
strumentation package at fixed known longitudinal 
positions within the collar, or by a free-falling in- 
strumentation package at known times or at known 
increments of time as the instrumentation package 
(^j moves through the collar. These measurements pro- 
^vide a longitudinal-position-dependent series of mag- 
netic field measurements BZ(z) which enable the 
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true magnitude of the terrestial magnetic field BZe in 
the direction of the longitudinal axis of the borehole 
to be calculated on the basis that BZ(z) = BZe + E- 
(z), where E(z) is the longitudinal-position-dependent 
longitudinal magnetic field error induced by the mag- 
netism of the drill string and the bottom-hole assem- 
bly. Several different methods of calculation are de- 
scribed, including polar and non-polar magnetic error 
function models. The methods can be extended to a 
full survey of the borehole heading by contempora- 
neous measurements of two further magnetic fields 
in each of two mutually orthogonal axes each also 
orthogonal to the longitudinal axis, along with con- 
temporaneous gravity vector component measure- 
ments in each of these three axes. Relevant meth- 
ods are described, along with apparatus for carrying 
out the heading survey methods. 
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Surveying 

This invention relates to the surveying of 
boreholes, and more particularly but not exclusively 
to determining the true azimuth of a borehole. 

When drilling a well for exploration and recov- 
ery of oil or gas, it is known to drill a deviated well, 
which is a well whose borehole intentionally de- 
parts from vertical by a significant extent over at 
least part of its depth. When a single drilling rig is 
offshore, a cluster of deviated wells drilled from 
that rig allows a wider area and a bigger volume to 
be tapped from the single drilling rig at one time 
and without expensive and time-consuming reloca- 
tion of the rig than by utilising only undeviated 
wells. Deviated wells also allow obstructions to be 
by-passed during drilling, by suitable control of the 
deviation of the borehole as it is drilled. However, 
to obtain the full potential benefits of well deviation 
requires precise knowledge of the instantaneous 
location and heading of the bottom-hole assembly 
(including the drilling bit and steering mechanisms 
such, as adjustable stabilisers). Depth of the 
bottom-hole assembly (or axial length of the 
borehole) can be determined from the surface, for 
example by counting the number of standard- 
length tubulars coupled into the drill string, or by 
less empirical procedures. However, determination 
of the location and heading of the bottom-hole 
assembly generally requires some form of down- 
hole measurement of heading. Integration of head- 
ing with respect to axial length of the borehole will 
give the borehole location relative to the drilling rig. 

In this context, the word "heading" is being 
used to denote the direction in which the bottom- 
hole assembly is pointing (ie. has its longitudinal 
axis aligned), both in a horizontal and vertical 
sense. Over any length of the borehole which can 
be considered as straight for the purposes of direc- 
tional analysis, the borehole axis in a deviated well 
will have a certain inclination with respect to true 
vertical. A vertical plane including this nominally 
straight length of borehole will have a certain angle 
(measured in a horizontal plane) with respect to a 
vertical plane including a standard direction; this 
standard direction is hereafter taken to be true 
magnetic north, and the said angle is the magnetic 
azimuth of the length of borehole under consider- 
ation (hereafter simply referred to as "azimuth"). 
The combination of inclination and azimuth at any 
point down the borehole is the heading of the 
borehole at that point; borehole heading can vary 
with depth as might be the case, for example, 
when drilling around an obstacle. 

Instrumentation packages are known, which 
can be incorporated in bottom-hole assemblies to 
measure gravity and magnetism in a number of 
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orthogonal directions related to the heading of the 
bottom-hole assembly. Mathematical manipulations 
of undistorted measurements of gravitational and 
magnetic vectors can produce results which are 

s representative of the true heading at the point at 
which the readings were taken. However, the mea- 
surements of magnetic vectors are susceptible to 
distortion, not least because of the masses of fer- 
rous materials incorporated in the drill string and 

70 bottom-hole assembly. Distortion of one or more 
magnetic vector measurements can give rise to 
unacceptable errors in the determination of head- 
ing, and undesirable consequences. Distortion of 
magnetic vectors in the region of the instrumenta- 

75 tion arising from inherent magnetism of conven- 
tional drill string and bottom-hole assembly compo- 
nents can be mitigated by locating the instrumenta- 
tion in a special section of drill string which is 
fabricated of non-magnetic alloy. However, such 

20 special non-magnetic drill string sections are rela- 
tively expensive. Moreover, the length of non-mag- 
netic section required to bring magnetic distortion 
down to an acceptable level increases significantly 
with increased mass of magnetic bottom-hole as- 

25 sembly and drill string components, with conse- 
quent high cost in wells which use such heavier 
equipment, eg. wells which are longer and/or deep- 
er. Hence such forms of passive error correction 
may be economically unacceptable. Active error 

30 correction by the mathematical manipulation of 
vector readings which are assumed to be error-free 
or to have errors which are small may give un- 
reliable results if the assumption is unwarranted. 
Before describing the invention, several defini- 

35 tions will be detailed with reference to Figs. 1 and 
2 of the accompanying drawings, wherein:- 

Fig. 1 is a schematic elevational view of the 
bottom-hole assembly of a drill string; and 

Fig. 2 is a schematic perspective view of 

40 various axes utilised for denoting directions in three 
dimensions. 

Referring first to Fig. 1 , the bottom-hole assem- 
bly of a drill string comprises a drilling bit 10 
coupled by a non-magnetic drill collar 1 2 and a set 
45 of drill collars 14 to a drill pipe 16. The drill collars 
14 may be fabricated of a magnetic material, but 
the drill collar 12 is substantially devoid of any self- 
magnetism. 

During local gravity and magnetic field vector 
50 measurements, the non-magnetic drill collar 12 
houses a downhole instrumentation package sche- 
matically depicted at 18. (In reality, the package 18 
would not be visible as is apparently the case in 
Fig. 1 since the package 18 is utilised within the 
interior of the collar 12). The downhole instrumen- 
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tation package 18 is capable of measuring gravity 
vectors and local magnetic vectors, for example by 
the use of accelerometers and fluxgates respec- 
tively. The instrumentation package 18 may be 
axially and rotationally fixed with respect to the 
bottom-hole assembly, including the drilling bit 10, 
whose heading is to be determined; the instrumen- 
tation package 18 would then be rigidly mounted in 
the bottom-hole assembly, within the non-magnetic 
drill collar 12 which is fabricated of non-magnetic 
alloy. Alternatively, the package 18 could be lower- 
ed through the collar 12, either on a wireline or as 
a free-falling package, with internal recording of the 
local gravity vectors and the local magnetic vec- 
tors. The alternative procedures for measurement 
processing according to whether the instrumenta- 
tion package 18 is axially fixed or mobile will be 
subsequently described. 

Referring now to Fig. 2 for convenience of 
conceptual presentation and calculation references, 
a hypothetical origin or omni-axial zero point "O" is 
deemed to exist in the centre of the instrumenta- 
tion package 18 (not shown in Fig. 2). Of the three 
orthogonal axes OX, OY and OZ defining the align- 
ment of the instrumentation relative to the bottom- 
hole assembly, the OZ axis lies along the axis of 
the bottom-hole assembly, in a direction towards 
the bottom of the assembly and the bottom of a 
borehole 20 drilled by the drilling bit 10. The OX 
and OY axes, which are orthogonal to the OZ axis 
and therefore lie in a plane 0.N2.E1 (now defined 
as the "Z-plane") at right angles to the bottom-hole 
assembly axis OZ, are fixed with respect to the 
body (including the collar 12) of the bottom-hole 
assembly. As viewed from above, the OX axis is 
the first of the fixed axes which lies clockwise of 
the upper edge of the (inclined) bottom-hole as- 
sembly, this upper edge lying in the true azimuth 
plane 0.N2.N1.V of the bottom-hole assembly. The 
angle N2.0.X in the Z-plane 0.N2.E1 (at right an- 
gles to OZ axis) between the bottom-hole assem- 
bly azimuth plane 0.N2.N1.V and the OX axis is the 
highside angle "HS". The OYaxis lies in the Z- 
plane 0.N2.E1 at right angles to the OX axis in a 
clockwise direction as viewed from above. A grav- 
ity vector measuring accelerometer (or other suit- 
able device) is fixedly aligned with each of the OX, 
OY and OZ axes. A magnetic vector measuring 
fluxgate (or other suitable device) is fixedly aligned 
in each of the OX, OY and OZ axes. The in- 
strumentation package 18 may be energised by 
any suitable known arrangement, and the instru- 
mentation readings may be telemetered directly or 
in coded form to a surface installation (normally the 
drilling rig) by any suitable known method, or alter- 
natively the instrumentation package 18 may incor- 
porate computation means to process instrumenta- 
tion readings and transmit computational results as 



distinct from raw data, or the instrumentation pack- 
age 18 may incorporate recording means for inter- 
nal recording of the local axial magnetic vectors for 
subsequent retrieval of the package 18 and on- 

s surface processing of the recorded measurements. 

Also.notionally vectored from the origin 0 are a 
true vertical (downwards) axis OV, a horizontal axis 
ON pointing horizontally to true Magnetic North, 
and an OE axis orthogonal to the OV and ON axes, 

ro the OE axis being at right angles clockwise in the 
horizontal plane as viewed from above (ie. the OE 
axis is a notional East-pointing axis). 

The vertical plane 0.N2.N1.V including the OZ 
axis and OV axis is the azimuth plane of the 

75 bottom-hole assembly. The angle V.O.Z between 
the OV axis and the OZ axis, ie. the angle in the 
bottom-hole assembly azimuth plane 0.N2.N1.V, is 
the bottom-hole assembly inclination angle "INC" 
which is the true deviation of the longitudinal axis 

20 of the bottom-hole assembly from vertical. Since 
the angles V.O.N1 and Z.O.N2 are both right an- 
gles and also lie in a common plane (the azimuth 
plane O.N2.N1.V), it follows that the angle N1.0.N2 
equals the angle V.O.Z, and hence "the angle 

25 N1 .0.N2 also equals the angle "INC". 

The vertical plane O.N.V. including the OV axis 
and the ON axis is the reference azimuth plane or 
true Magnetic North. The angle N.O.N1 measured 
in a horizontal plane O.N.N1.E.E1 between the ref- 

30 erence azimuth plane O.N.V. (including the OV axis 
and the ON axis) and the bottom-hole assembly 
azimuth plane O.N2.N1.V (including the OV axis 
and the OZ axis) is the bottom-hole assembly 
azimuth angle "AZ". 

35 The OX axis of the instrumentation package is 

related to the true Magnetic North axis ON by the 
vector sum of three angles as follows:- 

(1) horizontally from the ON axis round East- 
wards (clockwise as viewed from above) to* a hori- 

40 zontal axis O.N1 in the bottom-hole assembly azi- 
muth plane 0.N2.N1.V by the azimuth angle AZ 
(measured about the origin 0 in the horizontal , 
plane); 

* (2) vertically upwards from the horizontal 
45 axis O.N1 in the azimuth plane O.N2.N1.V to an 
inclined axis O.N 2 in the Z-plane (the inclined 
plane O.N2.E1 including the OX axis and the OY 
axis) by the inclination angle INC (measured about 
• the origin 0 in a vertical plane including the origin 

so 0); and 

(3) a further angle clockwise/Eastwards (as 
defined above) in the Z-plane from the azimuth 
plane to the OX axis by the- highside angle HS 
(measured about the origin 0 in the inclined Z- 
55 plane O.N2.E1 which includes the origin 0). 

Borehole surveying instruments measure the 
two traditional attitude angles, inclination and azi- 
muth, at points along the path of the borehole. The 
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inclination at such a point is the angle between the 
instrument longitudinal axis and the Earth's gravity 
vector direction (vertical) when the instrument lon- 
gitudinal axis as aligned with the borehole path at 
that point. Azimuth is the angle between the verti- 
cal plane which contains the instrument longitudinal 
axis and a vertical reference plane which may be 
either magnetically or gyroscopically defined; this 
invention is concerned with the measurement of 
azimuth defined by a vertical reference plane con- 
taining a defined magnetic field vector. 

Inclination and azimuth (magnetic) are conven- 
tionally determined from instruments which mea- 
sure the local gravity and magnetic field compo- 
nents along the directions of the orthogonal set of 
instrument-fixed axes [OX.OY.OZ]; traditionally. OZ 
is the instrument longitudinal axis. Thus, inclination 
and azimuth are determined as functions of the 
elements of the measurement set 
{GX,GY t GZ,BX,BY,BZ} where GX is the magnitude 
of the gravity vector component in direction OX.BX 
is the magnitude of the magnetic vector component 
in direction OX, etc. The calculations necessary to 
derive inclination and azimuth as .functions of 
GX,GY,GZ,BX t BY t BZ are well known. 

When the vertical magnetic reference plane is 
defined as containing the local magnetic field vec- 
tor at the instrument location, the corresponding 
azimuth angle is known as the raw azimuth; if the 
vertical magnetic reference plane is defined as 
containing the Earth's magnetic field vector at the 
instrument location, the corresponding azimuth an- 
gle is known as absolute azimuth. 

In practice, the value of the absolute azimuth is 
required and two methods to obtain it are presently 
employed: 

(i) The instrumentation package is contained 
within a non-magnetic drill collar (NMDC) which is 
sufficiently long to isolate the instrument from mag- 
netic effects caused by the proximity of the drill 
string (DS) above the instrument and the stabiliz- 
ers, bit, etc. forming the bottom-hole assembly 
(BHA) below the instrument In this case the 
Earth's magnetic field is uncorrupted by. the DS 
and BHA and the raw azimuth measured is equal to 
the absoluteazimuth. 

(ii) The corrupting magnetic effect of the DS 
and BHA is considered as an error vector along 
direction OZ thereby leaving BX and BY uncorrup- 
ted (components only of the Earth's magnetic 
field). The calculation of the absolute azimuth can 
then be performed as a function of 
GX.GY.GZ.BX.BY.Be, where Be is some value (or 
combination of values) associated with the Earth's 
magnetic field. 

The error in the measurement of absolute azi- 
muth by method_(ii) independent on the attitude of 



the measurement of the raw azimuth; the reasons 
for this are summarised as follows: 

(iii) the need to know the values of Earth's 
magnetic field components in instrument-magnetic- 

5 units to a high degree of accuracy: 

(iv) an inherent calculation error due to the 
availability of only the uncorrupted cross-axis 
(BOXY) magnetic vector component. [This is an- 
alagous to measuring only the gravity component 

w GZ and then attempting to determine the inclination 
(INC) from INC = ACOS(GZ), with the magnitude 
of Earth's gravity = 1 instrument gravity-unit]. 

It is therefore an object of the invention to 
provide an improved method of surveying a 

75 borehole, and more particularly but not exclusively 
to provide an improved method of surveying the 
magnetic azimuth of a borehole. 

According to the present invention there is 
provided a method of surveying the magnetic azi- 

20 muth of a borehole penetrated by a magnetic drill 
string coupled through a substantially non-magnetic 
drill collar to a magnetic bottom-hole assembly, by 
deriving the true magnitude of the terrestial mag- 
netic field BZe in the direction of the longitudinal 

25 axis of the borehole (the OZ axis as defined) in the 
region of the substantially non-magnetic drill collar, 
said method comprising the steps of measuring the 
longitudinal magnetic field BZ (the component of 
the magnetic field B in the direction OZ) at a 

30 plurality of points along the length of the substan- 
tially non-magnetic drill collar to provide a 
longitudinal-position-dependent series of magnetic 
field measurements BZ(z), and calculating BZe on 
the basis that BZ(z) = BZe + E(z), where E(z) is 

35 the longitudinal-position-dependent longitudinal 
magnetic field error induced by magnetism of the 
drill string and bottom-hole assembly. 

The calculation of BZe may be based on the 
assumption that the longitudinal magnetic field er- 

40 ror E(z) is induced by a plurality of notional mag- 
netic poles longitudinally distributed along the lon- 
gitudinal axis adjacent the substantially non-mag- 
netic drill collar. The plurality of notional magnetic 
poles assumed to be inducing the longitudinal 

45 magnetic field error E(z) may comprise one pole 
pair or a plurality of pole pairs. 

However, it is not essential to calculate the 
longitudinal magnetic field error E(z) in terms of a 
magnetic pole model; any mathematical method or 

so curve-matching exercise which results in the gen- 
eration of a function E(z) such that the measured 
distribution BZ(z) is closely represented by E(z) + 
K (where K is a constant) is sufficient to determine 
BZe = K. 

55 Moreover, in order to generate the longitudinal- 
position-dependent longitudinal magnetic error E- 
(z). it is not necessary to know the absolute posi- 
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surements of longitudinal magnetic field BZ are 
made to provide the longitudinal-position-depen- 
dent series of magnetic field measurement BZ(z). It 
is sufficient to know the positions of the measure- 
ment points relative to each other in order to deter- 5 
mine the longitudinal magnetic error E(z). 

The relative positions of the measurement 
points are known for the case where the instrumen- 
tation package or other local axial magnetic field 
vector measuring means contains a plurality of OZ w 
fluxgates at known mutual spacings along the lon- 
gitudinal 2 axis and is static within the NMDC at 
the time of measurement, and also for the case 
where the instrumentation package or other mea- 
suring means is suspended from a wireline and is 
passes longitudinally through the non-magnetic drill 
collar at known depths controlled from the surface 
above the well. 

In the case where the instrumentation package 
falls freely through the non-magnetic drill collar, 20 
measurements are generally not made at known 
increments of distance, but are made (and re- 
corded) at known times or at known increments of 
time; a procedure for converting such time-sepa- 
rated measurements to distance-separated mea- 25 
surements is also comprised within the scope of 
the present invention and will be described subse- 
quently. 

The present invention also provides apparatus 
for carrying out the foregoing magnetic azimuth 30 
surveying method, said apparatus comprising an 
instrumentation package containing at least two 
longitudinal magnetic field measuring devices hav- 
ing a known fixed separation(s). 

Said apparatus may alternatively comprise an 35 
instrumentation package containing at least two 
longitudinal magnetic field measuring devices hav- 
ing a known fixed mutual separation(s), and a re- 
cording means to which said magnetic field mea- 
suring devices are connected for recording a plu- *o 
rality of longitudinal magnetic field measurements 
performed by each said device at known times or 
at known increments of time as said instrumenta- 
tion package moves through said substantially non- 
magnetic drill collar. 

The foregoing magnetic azimuth surveying 
method may be extended to provide a method of 
surveying the heading of a borehole penetrated by 
a magnetic drill string coupled through a substan- 
tially non-magnetic drill collar to a magnetic so 
bottom-hole assembly, by deriving the true mag- 
nitude of the terrestial magnetic field BZe in the 
direction of the longitudinal axis of the borehole in 
the region of the substantially non-magnetic drill 
collar, said method comprising the steps of mea- 55 
suring the longitudinal magnetic field BZ at a plu- 
rality of points along the length of the substantially 
non-magnetic drill collar to provide a longitudinal- 



position-dependent series of magnetic field mea- 
surements BZ(z), contemporaneously measuring 
the magnetic fields Bx and By in two mutually 
orthogonal axes each also orthogonal to the lon- 
gitudinal axis, contemporaneously measuring grav- 
ity vector components in each of the said three 
axes to produce respective gravity vector measure- 
ments Gx, Gy and Gz, calculating BZe on the basis 
that BZ(z) = BZe + E(z), where E(z) is the 
longitudinal-position-dependent longitudinal mag- 
netic field error induced by magnetism of the drill 
string and the bottom-hole assembly, and solving 
the function [Gx, Gy, Gz, Bx, By. BZe] to deter- 
mine said heading. 

The present invention further provides appara- 
tus for carrying out the immediately foregoing 
method of surveying the heading of a borehole, 
said apparatus comprising an instrumentation pack- 
age containing at least two longitudinal magnetic 
field measuring devices having a known fixed mu- 
tual separation(s), two further magnetic field mea- 
suring devices for contemporaneously measuring 
magnetic fields in two mutually orthogonal axes 
each also orthogonal to the longitudinal axis, three 
gravity vector component measuring devices for 
contemporaneously measuring gravity vector com- 
ponents in each of the said three axes, and a 
recording means to which each of said magnetic 
field measuring devices and each of said gravity 
vector component measuring devices is connected 
for recording the respective measurements of the 
respective magnetic fields and the respective mea- 
surements of the respective gravity vector compo- 
nents when said intrumentation package is within 
the substantially non-magnetic drill collar. 

The present invention still further provides ap- 
paratus for carrying out the method of surveying 
magnetic azimuth and for carrying out at least the 
magnetic azimuth survey step of the method of 
surveying the heading, of a borehole, said appara- 
tus comprising an instrumentation package contain- 
ing longitudinal magnetic measuring means for 
measuring the longitudinal magnetic field at a plu- 
rality of positions along the longitudinal axis, said 
instrumentation package further containing deter- 
mining means for directly or indirectly determining 
the respective absolute or relative distances along 
the longitudinal axis of the positions at which the 
plurality of longitudinal magnetic field measure- 
ments are made. 

Embodiments of the invention will now' be de- 
scribed by way of example with reference to Figs. 
3-9 of the accompanying drawings wherein: 

Fig. 3 is a graphical representation of the 
variation of azimuth reading errors with inclination, 
for a typical present day instrumentation package; 

Fig. 4 is a schematic representation of a 
simple model of an error-inducing notional mag- 
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netic pole system; 

Fig. 5 is a schematic representation of a 
complex model of an error-inducing notional mag- 
netic pole system; 

Fig. 6 is a graphical representation of cal- 
culated results employing one model of field sys- 
tem; 

Fig. 7 is a graphical representation of cal- 
culated results employing another model of field 
system; 

Fig. 8 is a schematic representation of a 
free-fall instrumentation package for measuring and 
recording local longitudinal magnetic fields at 
points having a fixed known mutual separation and 
at known times or at known increments of time; 
and 

Fig. 9 is a graphical representation of part of 
a procedure for converting the time-separated mea- 
surement obtained by the instrumentation of Fig. 8 
to distance-separated measurements. 

Fig. 3 indicates the relative accuracies of deter- 
mining the raw and absolute azimuths for the 
worst-case situation when the local axial magnetic 
field vector measuring instrument is lying with its 
longitudinal axis east/west; the values are calcu- 
lated using a set of errors representative of the 
limit of what is achievable for present-day instru- 
ments. 

For the sake of convenience in referring to 
magnetic perturbations of local magnetic fields, the 
bottom hole assembly comprising the drilling bit 10 
(and any associated magnetic components) will 
subsequently be referred to as the "BHA", the drill 
collars 14 and drill string 16 (plus any associated 
magnetic components) will subsequently be re- 
ferred to as the "DS", and the non-magnetic drill 
collar 12 will subsequently be referred to as the 
"NMDC". 

This invention concerns a method of determin- 
ing absolute azimuth without the need to use ac- 
curate Earth's field data and without the problems 
associated with degradation of the calculation due 
to attitude changes. The method itself is dependent 
on two key factors: 

(a) Acknowledge of the axial magnetic com- 
ponent (BZ) at distributed points along the axis of 
the NMDC. 

(b) The selection of a theoretical magnetic 
model to represent the cause of the corrupting field 
due to effects of the DS and BHA. 

The accuracy of the method is entirely depen- 
dent on the extent to which data of (a) is known 
since this determines the degree of sophistication 
which can be used to select the model (b). (While 
the method will first be described in terms of a 
magnetic pole model, variations of the method em- 



Unless the length of the NMDC is very small, 
say less than 10 feet, experience shows (as might 
be expected from Fig. 1) that the effect of the 
magnetic DS and BHA material is to produce a 

s magnetic error field {E} at a point on the axis of 
the NMDC and remote from its end with the direc- 
tion of (E) substantially along the longitudinal axis 
OZ of the NMDC. 

Magnetic models representing this magnetic 

io configuration can be reasonably postulated in 
terms of notional magnetic poles of various 
strengths distributed along the NMDC axis (OZ) 
direction. The degree of sophistication for such 
models will be dependent upon both the number of 

is such magnetic poles employed and the degrees of 
freedom in their positioning. 

The principles of the magnetic polar models 
will now be described. 

If the value of BZ is measured at various points 

20 z as measured along the NMDC length (z being 
zero at one end of the NMDC), then at any such 
point z, the longitudinal-position-dependent value of 
BZ is BZ(z) such that:- 
BZ(z) = BZe + E(z) 

25 where BZe is the value of the Earth's magnetic 
field component along OZ t and E(z) is the 
longitudinal-position-dependent value of the error 
field {E} at that point. 

In terms of any postulated polar magnetic 

30 model, E(z) will be a function both of the notional 
pole strengths and of distances (functions of z) 
from the points of the notional poles employed in 
the model, but BZe is invariant with respect to z. If 
measurement of BZ(z) are made at points along 

35 the OZ axis inside the NMDC, then sets of equa- 
tions can be formed and solved for the unknowns 
of the model as well as for BZe. Clearly the num- 
ber of unknowns for the model which can be deter- 
mined in this manner will be dependent on the 

40 number of equations so formed; ie. on the number 
of points at which BZ(z) is measured along the 
NMDC length. 

Some magnetic polar models will now be de- 
scribed in detail. 

45 Examples of two magnetic polar models are 
considered here; the first example is the simplest 
possible configuration of magnetic poles which 
might be employed and the second example is 
probably beyond the limit to which the sophistica- 

50 tion for such models needs to be taken to produce 
more accurate results. 



(i) The simple model is schematically depicted in 
55 Fig. 4. 

r^rc -n,Hp' rn*T* i r*r»c l ♦Hq* the effect of the DS 
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pole strengths located at each end of the NMDC. 
each pole having a longitudinal field strength P. 

The value of the axial field at distance z from 
the upper end of the NMDC can be written in terms 
of this model as: 

BZ(2) = BZe + E(z) = BZe + P/fl-z) 2 
with unknowns BZe and P. 

Clearly, if measurements are made at two 
points along the NMDC axis, then two such equa- 
tions are obtained which can be used to solve for 
BZe and P (in instrument-magnetic-units). It should 
be noted that the selection of the locations of the 
two points at which the measurements are made 
will be important in practice. 



falling and be retrieved when the complete string is 
pulled from the hole, or. alternatively, a wireline 
may be used both to lower and to retrieve the SIA]. 

With present-day survey instruments, measure- 
ments of BZ(z) could be made at short time inter- 
vals and stored in memory. The data recorded as 
the SIA leaves the DS and transverses the NMDC 
can be correlated with distance along the NMDC 
axis for a known or presumed velocity profile or 
constant velocity and. thus, the BZ(z) profile for 
this transverse can be stored for future processing 
to determine the magnetic pole model characteris- 
tics necessary to allow the determination of BZe at 
the SIA location. 



(ii) The complex model-is schematically depicted in 
Fig. 5. 

This model considers the. effect of the DS and 
BHA in terms of four poles with pole strengths 
P1.P2.P3 and P4 located at distances L1.L2.L3 and 
L4 respectively from the upper end of the NMDC. 

The value of the axial field at distance z from 
the upper end of the NMDC can be written in terms 
of this model as: 
BZ(z) = BZe + E(z) 

where _ 
E(z) = - P1/(L1+z) 2 + P2/(L2 + z) 2 + P3/(L3-z) 2 - 

P4/(L4-z) 2 

The unknowns in this case are 
P1.P2.P3,P4.L1,L2,L3.L4 and Bze. Clearly, at least 
9 measurements of BZ(z) must be made in order 
to fully characterise this model. 



Acquisition of Data 

In order to determine the characteristics of the 
Error Function E(z) generation model used to pre- 
dict the effects of the DS and BHA on the magnetic 
field at points within the NMDC, it is necessary to 
measure the total axial magnetic field component 
BZ(z) at points along the axis of the NMDC. Clear- 
ly, the-instrument package could consist of a series 
of axial fluxgates at appropriate spacings in addi- 
tion to the normal configuration of three gravity 
sensors plus three magnetic fluxgates. However, 
there are ways to obtain the BZ(z) profile points 
without the necessity to change to any great extent 
the present surveying operational procedures: 



(i) Single-Shot Survey:- 

A survey instrument assembly (SIA) is passed 
,^...„ thrn..nh the ns tn a known location within the 



(ii) Multishot survey:- 

The SIA, which normally contains at least two 

20 magnetic survey instruments, is free-dropped to a 
known location in the NMDC. Again. BZ(z) can be 
measured and stored as the SIA transverses the 
NMDC to its location; with the multiplicity of survey 
instrument data, it is possible to characterise ac- 

25 curately an Error Function E(z) generation model 
representative of the DS and BHA at the bottom- 
hole location. 

Survey instrument(s) data is then recorded as 
the complete string assembly is pulled from the 

ao hole- it is possible that, due to induced magnetisa- 
tion effects, the parameters of the model will need 
revision as the attitude of the NMDC and SIA 
changers. For examples at any survey point, the 
pole strengths in a magnetic pole model can be 

os scaled according to the difference in BZ from two 
survey instruments spaced at appropriate points 
along the NMDC axis. Thus, using these models, 
BZe values can be determined for each survey 

point. ... . 

40 (A procedure for determining BZe by utilising 
two axial fluxgates or equivalent devices in a free- 
fall SIA, and which is applicable to polar and non- 
polar models, is detailed subsequently). 



45 



Relative Accuracies 



In the discussion that follows 1 instrument- 
magnetic-unit is approximately equal to 1 micro- 

so tesla. The determination of the Earth's magnetic 
field component BZe in instrument-magnetic-units 
from the Error Function E(z) generation model is 
dependent on the degree to which the model cho- 
sen is representative of the DS and BHA effects 

55 and the accuracy to which differences in BZ(z) at 
points along the axis (OZ) of the NMDC can be 
measured; with a multiplicity of data points along 
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model with sufficient sophistication to represent 
very closely DS and BHA magnetic effects, and 
differences in BZ(z) values along the NMDC axis 
will be independent of the 02-fluxgate datum er- 
rors. Therefore, since it should be possible in prac- 
tice to match fluxgate scale factors for the 
OX.OY.OZ fluxgates within an error band of width 
+/-0.1%, it is reasonable to suppose that the error 
band for BZe derived from model approach could 
be better than +/- 0.2 instrument-magnetic-units. 

Methods which derive absolute azimuth as a 
function of (GX,GY,GZ,BX,BY ( Be'), where Be' is an 
assumed known value of one or more components 
of the Earth's magnetic field Be at the drilling 
location, effectively require measurements of the 
magnetic field components (BX.BY.Be) in absolute 
units. Given this necessity to match the scale fac- 
tors of the survey instrument fluxgates to an ab- 
solute reference, it is optimistic to assume that any 
component value of the Earth's magnetic field used 
in the calculation can be known in practice to an 
accuracy of better than + /- 0.2 instrument- 
magnetic-units. 

Fig. 6 shows a comparison for the Error Func- 
tion E(z) generation model method of this invention 
and a calculation which determines absolute azi- 
muth as a function of (GX.GY.GZ.BY.BVe), where 
BVe is the value of the vertical component of the 
Earth's magnetic field Be at the drilling location 
(assumed known from independent sources). The 
error in BVe is taken as + /-0.2 instrument- 
magnetic-units (optimistic) and the error in BZe 
from the model method is taken as 0.4 instrument- 
magnetic-units (pessimistic). The value of the ab- 
solute (or raw) azimuth which would be obtained in 
a long NMDC configuration with the same instru- 
ment error set is also plotted. 

The results are based on instrument error' 
bands as follows: 

Gravity sensors: Scale factors +/-0.1% 
Datums +/-0.1%*g 

Magnetic sensors: Scale factors +/-0.2% 
Datums +/-0.2%'Be 

100 sample calculations are performed for 
each inclination value with the true azimuth taken 
as 90 degrees (east); the instrument error set and 
the instrument rotation angle (about OZ) are ran- 
domly chosen for each calculation. For comparison 
purposes, the absolute value of the mean error plus 
twice the standard deviation is the parameter plot- 
ted. 

Fig. 7 shows the same plots for calculations 
with the magnetic sensor's scale factor error re- 
duced to +/-0.10% and the error in BZe from the 
model method taken as 0.2 instrument-magnetic- 
units. 



Concluding Comments On Polar Models:- 

With present-day survey .instruments capable 
of measuring and recording the BZ(z) component 

5 of the local magnetic field within the NMDC at a 
frequency of several times per second, it is possi- 
ble to obtain a highly detailed profile of the axial 
magnetic field within the NMDC. The profile can be 
used to characterise an axial magnetic pole dis- 

io tribution model which will represent the magnetic 
effect of DS and BHA at points along the axis of 
the NMDC to a high degree of accuracy. Using this 
model, the corrupting field can be estimated at any 
point along the axis of the NMDC and. thus, the 

75 axial component of the Earth's magnetic field (BZe) 
can also be estimated at any such point. 

The results of calculations performed and sum- 
marised in the plots of Figs. 6 and 7 suggest that 
this method will be much superior to the currently 

20 used calculations which require an accurate knowl- 
edge of the Earth's magnetic field from an in- 
dependent source. While there is probably little to 
choose between the methods at inclinations up to 
about 40 degrees, at greater inclinations the polar 

25 method is likely to yield much better results. 

Clearly, the most accurate method of obtaining 
absolute azimuth is still through the employment of 
a (sufficiently) long NMDC to minimize the DS and 
BHA effects, but length and cost considerations do 

30 not necessarily make this the most attractive 
means of measurement and the operational advan- 
tages of running with a shorter NMDC are consid- 
erable. 

35 

Non-polar Derivations of BZe:- 

lt has been described above how measurement 
of BZ(z) can be made at a sufficient number of 
40 points along the NMDC axis to permit the solution 
of a set of simultaneous equations, each in the 
form :- 

BZ(z) = BZe + E(z) 

such as the yield the OZ vector value BZe of the 
45 Earth's magnetic field (which is the objective of the 
procedure). The minimum number of such mea- 
surements is determined by the complexity of the 
magnetic pole model used to generate the mag- 
netic distortion function E(z). 
50 However, it is not essential to calculate the 

function E(z) in terms of a magnetic pole model; 
any mathematical method which results in the gen- 
eration of a function E(z) such that the measured 
distribution BZ(z) is closely represented by E(z) + 
55 K (where K is a constant) is sufficient to determine 
BZe = K. 

Moreover, in order to aenerate the longitudinal- 
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(z), it is not necessary to know the absolute posi- 
tions of the measurement points at which the mea- 
surements of longitudinal magnetic field BZ are 
made to provide the longitudinal-position-depen- 
dent series of magnetic field measurements BZ(z). 
It is sufficient to know the positions of the measure- 
ment points relative to each other in order to deter- 
mine the longitudinal magnetic error E(z). 

The relative ppsitions of the measurement 
points are known for the case where the instrumen- 
tation package 18 or other local axial magnetic field 
vector measuring survey-instrument assembly SIA 
contains a plurality of OZ fluxgates (or of equiv- 
alent magnetic measuring devices) at known mu- 
tual spacings along the longitudinal OZ axis and is 
static within the NMDC 12 at the time of measure- 
ment, and also for the case where the instrumenta- 
tion package 18 or other SIA is suspended from a 
wireline and passes longitudinally through the 
NMDC 12 at a velocity controlled by the wireline 
operator on the surface above the well. 

In the case where the instrumentation package 
18 or other SIA falls freely through the NMDC 12, 
measurements are not made at known increments 
of distance because of the uncontrolled rate of fall, 
but are made and recorded at known times or at 
known increments of time. A modified form of 
intrumentation package 18 and a procedure of con- 
verting such time-separated measurements for sub- 
sequent calculation of BZe will now be described 
with reference to Figs. 8 and 9. 

Referring first to Fig. 8, the modified instru- 
mentation package 18 comprises a first local axial 
(OZ) vector measuring fluxgate F1 mounted at the 
upper (trailing) end of the package 18, and a sec- 
ond local axial (OZ) vector measuring fluxgate F2 
mounted at the lower (leading) end of the package 
18. The fluxgates F1 and F2 have a fixed mutual 
axial separation 5d ('delta-d') within the package 
18. Both fluxgates F1 and F2 are connected to an 
internal recording device Rec. which records fre- 
quent B(z) measurements at known increments of 
time (or in any other time-dependent reproducible 
manner) as the package 18 free-falls through the 
NMDC A2. The instrumentation package 18 also 
includes fluxgates Fx and Fy respectively measur- 
ing the local magnetic field vectors in the OX and 
OY directions, as well as local gravity vector mea- 
suring accelerometers Gx,Gy and Gz, respectively 
for measuring the local gravitational vector Gx 
along the OX axis, for measuring the local gravita- 
tional vector Gy along the OY axis, and for measur- 
ing the local gravitational vector Gz along the OX 
axis. The fluxgates Fx and Fy, and the accelero- 
meters Gx.Gy and Gz are also connected to the 
internal recording device Rec. so as to make local 
gravity vector measurements correlated in time, 
and hence in - position, with the local magnetic 



vector measurements. 

Referring now to Fig. 9. this shows a twin 
graph of the two plots of the time-dependent local 
longitudinal (OZ) magnetic vector BZ(t) with re- 

s spect to time T as measured by each , of the 
fluxgates F1 and F2 (and recorded in the recorder 
Rec.) while the instrumentation package 18 freely 
falls down through the NMDC 12. Individual record- 
ings are not denoted on either plot, the discrete 

w markings being subsequently added at selected 
pairs of points,, one on each plot, which are of 
mutually equal values of BZ(t), though not nec- 
essarily at any particular values of BZ(t). The rea- 
sons for the addition of such markings are given 

15 below. 

Taking either of the individual plots of BZ(t) in 
Fig. 9, the valley-shaped plot is characteristic of 
the longitudinal magnetic field vector diminishing 
from an initially high value of BZ as the respective 

20 fluxgate leaves the drill string DS and its immediate 
local magnetic influence, falling to a non-zero mini- 
mum approximately mid-way between the drill 
string DS and the bottom-hole assembly BHA, and 
rising again as the instantaneous BZ is increasingly 

25 influenced by the approach of the fluxgate to the 
BHA with its local magnetic influence. If the in- 
strumentation package 18 falls at a substantially 
constant velocity, the two plots will be substantially 
identical, but mutually slightly displaced along the 

30 horizontal time axis V, whereas if the package 18 
changes its velocity [due to transient or continuous 
acceleration(s) and/or deceleration^)], the two plots 
will not be identical. However, the procedure de- 
scribed below enables the time-dependent plots 

35 BZ(t) to be converted to the requisite position- 
dependent plots BZ(z) for subsequent calculation 
of BZe, without any need to assume any particular 
constant velocity or velocity profile for the instru- 
mentation package 18 in its uncontrolled lorigitudi- 

40 nal passage through the NMDC 12. (The procedure 
is also applicable to the case where the instrumen- 
tation package 18 is lowered at a known or con- 
trolled velocity (eg. by being lowered on a wireline) 
but such known or controlled velocity does not 

45 have to be taken into account). 

The time/position conversion procedure de- 
pends on the fact that, regardless of velocity or of. 
velocity changes, each of the fluxgates Ft and F2 
will pass through the same longitudinal position 

so along the OZ axis (albeit at different times), and 
hence through the same local longitudinal magnetic 
field. Thus any two adjacent points, one on each 
adjacent fluxgate plot, which are at mutually equal 
values of the local longitudinal magnetic field BZ(t) 

55 represent the successive passages of the two flux- 
gates through the same longitudinal position. The 
horizontal separations of any such adjacent pair of 
equi-valued points of BZ(t) is the time interval 5t- 
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( 'delta-t') from the passage of the leading fluxgate 
F2 until the trailing fluxgate F1 passes the same 
point. Since the fluxgates F1 and F2 have a known 
separation 66 which is constant (invariant with re- 
spect to time), this separation 5d divided by the 
relevant time interval 6t at any point of the traversal 
of the NMDC 12 is the velocity of the package 18 
at that point. This yields a velocity/time profile 
which can be integrated to derive distance values 
giving relative position's at which the intially se- 
lected values of BZ apply. 

Reverting to Fig. 9 t adjacent pairs of points on 
the two plots of BZ(t) are selected, at mutually 
identical values of BZ(t). The points on the plot of 
measurements form the trailing fluxgate F1 are 
denoted by a " + ", while the points on the plot of 
measurements from the leading, fluxgate F2 are 
denoted by an "o". For any arbitarily selected point 
on one plot, there is a unique adjacent point on the 
adjacent plot at the identical value of BZ(t). me 
actually selected points need not have any specific 
value, nor any mutually related values, save that 
their number and distribution are at least sufficient 
to provide the requisite accuracy in producing the 
resultant velocity/time profile; by way of example 
only, Fig. 9 depicts eight such pairs of points at 
approximately equal intervals along the horizontal 
time axis t . 

The time Y attributed to any given pair of 
points on the pair of BZ(t) curves can be referen- 
ced to the trailing fluxgate F1 (points denoted 
" + "), or referenced to the leading fluxgate F2 
(points denoted "o"), or referenced to a point mid- 
way between these points, as illustrated by way of 
example in Fig. 9 for the second pair of points 
only. 

Having obtained the speed/time function and 
then (by integration) obtained the distance/time 
function therefrom, as the basis of derivable rela- 
tive positions, all as described above, the resultant 
derived values of BZ(z) can be utilised in any 
suitable polar or non-polar magnetic error function 
model as previously described to derive the value 
of BZe as the value of the longitudinal (OZ axis) 
vector component of the terrestial. magnetic field 
within the borehole 20 at the time and place of the 
original measurements of local magnetic and grav- 
ity vectors. This value of BZe, in conjunction with 
the comtemporaneous measured values 
Bx,By,Gx,Gy and Gz of the local gravity vectors 
(produced respectively by the fluxgates Fx and Fy, 
and the accelerometers Gx.Gy and Gz within the 
modified instrumentation package 18 of Fig. 8), 
yield a function (Gx.Gy.Gz.Bx.By.Bz) which can be 
resolved as previously described to yield the head- 
ing of the borehole 20 at the location of the NMDC 
12. 

Alternative orocedyres for resolving the 



borehole heading from fewer than all six orthogonal 
gravity and magnetic vectors may be employed 
without departing from the scope of the invention, 
which essentially lies in/the novel method of deter- 

5 mining BZe. It is equally within the scope of the 
. present invention that if the value of BZe were the 
only unknown to be determined, this single un- 
known could be determined by the method of the 
present invention. (In either of these cases, one or 

w more of the fluxgates Fx and Fy and/or the ac- 
celerometers Gx,Gy and Gz within the instrumenta- 
tion package 18 of Fig. 8 might then be redundant, 
but this would not affect the essential scope of the 
method of the present invention). 

75 While certain modifications and variations have 
been described above, the invention is not re- 
stricted thereto, and other modifications or vari- 
ations can be adopted without departing from the 
scope of the invention as defined in the appended 

20 Claims. 



Claims 

25 1 . A method of surveying the magnetic azimuth 

of a borehole penetrated by a . magnetic drill string 
coupled through a substantially non-magnetic drill 
collar to a magnetic bottom-hole assembly, by de- 
riving the true magnitude of the terrestial magnetic 

30 field BZe in the direction of the longitudinal axis of 
the borehole in the region of the substantially non- 
magnetic drill collar, said method comprising the 
steps of measuring the longitudinal magnetic field 
BZ at a plurality of points along the length of the 

35 substantially non-magnetic drill collar to provide a 
longitudinal-position-dependent series of magnetic 
field measurements BZ(z), and calculating BZe on 
the basis that BZ(z) = BZe + E(z), where E(z) is 
the longitudinal-position-dependent longitudinal 

40 magnetic field error induced by magnetism of the 
drill string and the bottom-hole assembly. 

2. A method as claimed in Claim 1 wherein the 
calculation of BZe is based on the assumption that 
the longitudinal magnetic field error E(z) is induced 

45 by a plurality of notional magnetic poles longitudi- 
nally distributed along the longitudinal axis adjacent 
the substantially non-magnetic drill collar. 

3. A method as claimed in Claim 2 wherein the 
plurality of notional magnetic poles assumed to be 

so inducing the longitudinal magnetic field error E(z) 
comprises one pole pair or a plurality of pole pairs. 

4. A method as claimed in Claim 1 comprising 
the step of generating a function E(z) such that the 
measured distribution BZ(z) is closely represented 

55 by E(z) + K (where K is a constant) such that BZe 
= K. 

5. A method as claimed in any of Claims 1 to 5 
wherein the measurements of longitudinal magnetic 
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field BZ are performed by at least two longitudinal 
magnetic field measuring devices having a known 
fixed mutual separation(s) and which are passing 
longitudinally through the substantially non-mag- 
netic drill collar during said measurements, a plu- 
rality of such measurements being performed by 
each said device at known times or at known 
increments of time to produce respective time- 
dependent local longitudinal magnetic field vectors 
BZ(t), deriving increments of time therefrom at se- 
lected values of BZ on the basis that said devices 
successively pass through any given longitudinal 
position and measure equal values of BZ thereat 
such that said increments of time represent the 
time differences of such successive passes, divid- 
ing said increments of time by the mutual 
separation(s) of said, devices to derive a 
velocity/time function of the passage of said de- 
vices through said substantially non-magnetic drill 
collar, and integrating said velocity/time function to 
derive distance values giving relative positions at 
which the selected values of BZ apply whereby to 
derive said longitudinal-position-dependent series 
of magnetic field measurements BZ(z). 

6. A method of surveying the heading of a 
borehole penetrated by a magnetic drill string 
coupled through a substantially non-magnetic drill 
collar to a magnetic bottom-hole assembly, by de- 
riving the true magnitude of the terrestial magnetic 
field BZe in the direction of the longitudinal axis of 
the borehole in the region of the substantially non- 
magnetic drill collar, said method comprising the 
steps of measuring the longitudinal magnetic field 
BZ at a plurality of points along the length of the 
substantially non-magnetic drill collar to provide a 
longitudinal-position-dependent series of magnetic 
field measurements BZ(z), contemporaneously 
measuring the magnetic fields Bx and By in two 
mutually orthogonal axes each also orthogonal to 
the longitudinal axis, contemporaneously measur- 
ing gravity vector components in each of the said 
three axes to produce respective gravity vector 
measurements Gx. Gy and Gz, calculating BZe on 
the basis that BZ(z) = BZe + E(z), where E(z) is 
the lojpgitudinal-position-dependent longitudinal 
magnetic field error induced by magnetism of the 
drill string and the bottom-hole assembly, and solv- 
ing the function [Gx, Gy, Gz, Bx, By, BZe] to 
determine said heading. 

7. Apparatus for carrying out the method as 
claimed in any of Claims 1 to 5, said apparatus 
comprising an instrumentation package containing 
at least two longitudinal magnetic field measuring 
devices having a known fixed mutual separation(s). 

8. Apparatus for carrying out the method as 
claimed in Claim 5, said apparatus comprising an 



ing a known fixed mutual separation(s), and a re- 
cording means to which said magnetic field mea- 
suring devices are connected for recording a plu- 
rality of longitudinal magnetic field measurements 
s performed by each said device at known times or 
at known increments of time as said instrumenta- 
tion package moves through said substantially non- 
magnetic drill collar. 

9. Apparatus for carrying out the method as 
w claimed in Claim 6, said apparatus comprising an 

instrumentation package containing at least two 
longitudinal magnetic field measuring devices hav- 
ing a known fixed mutual separation (s), two further 
magnetic field measuring devices for contempora- 

75 neousiy measuring magnetic fields in two mutually 
orthogonal axes each also orthogonal to the lon- 
gitudinal axis, three gravity vector component mea- 
suring devices for contemporaneously measuring 
gravity vector components in each of the said three 

20 axes, and a recording means to which each of said 
magnetic field measuring devices and each of said 
gravity vector component measuring devices is 
connected for recording the respective measure- 
ments of the respective magnetic fields and the 

25 respective measurements of the respective gravity 
vector components when said intrumentation pack- 
age is within the substantially non-magnetic drill 
collar. 

10. Apparatus for carrying out the method as 
30 claimed in any of Claims 1 to 6, said apparatus 

comprising an instrumentation package containing 
longitudinal magnetic measuring means for mea- 
suring the longitudinal magnetic field at a plurality 
of positions along the longitudinal axis, said in- 

35 strumentation package further containing determin- 
ing means for directly or indirectly determining the 
respective absolute or relative distances along the 
longitudinal axis of the positions at which the plu- 
rality of longitudinal magnetic field measurements 

40 are made. 
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0 Borehole surveying methods and apparatus for 
surveying the true longitudinal magnetic field within 
a substantially non-magnetic drill collar occupying 
the part of a borehole being surveyed, despite the 
collar being of insufficient length to provide longitudi- 
nal magnetic field measurements which are uncor- 
rupted by the longitudinal magnetic influences of 
adjacent magnetic drill string and bottom-hole as- 
sembly components. A plurality of longitudinal mag- 
netic field measurements are made by a static in- 
strumentation package at fixed known longitudinal 
positions within the collar, or by a free-falling in- 
strumentation package at known times or at known 
increments of time as the instrumentation package 
moves through the collar. These measurements pro- 
vide a longitudinal-position-dependent series of mag- 
netic field measurements BZ(z) which enable the 



true magnitude of the terrestial magnetic field BZe in 
the direction of the longitudinal axis of the borehole 
to be calculated on the basis that BZ(z) = BZe + E- 
(z), where E(z) is the longitudinal-position-dependent 
longitudinal magnetic field error induced by the mag- 
netism of the drill string and the bottom-hole assem- 
bly. Several different methods of calculation are de- 
scribed, including polar and non-polar magnetic error 
function models. The methods can be extended to a 
full survey of the borehole, heading by contempora- 
neous measurements of two further magnetic fields 
in each of two mutually orthogonal axes each also 
orthogonal to the longitudinal axis, along with con- 
temporaneous gravity vector component measure- 
ments in each of these three axes. Relevant meth- 
ods are described, along with apparatus for carrying 
out the heading survey methods. 
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